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ABSTRACT 


For the past eight years, the small automated telescope Raven has been 
tested in detecting and tracking deep space objects. As the Raven has proven 
successful in tracking this regular and predictable orbit, its one arc-second 
accuracy made it a perfect candidate to attempt to accurately track the less 
predictable Highly Eccentric Orbit (HEO) objects. Ranging data was obtained 
from the Sirius satellite radio company for the Sirius3 satellite (Satellite Control 
Center (SCC) # 26626). This satellite was chosen for its long dwell time over the 
United States and for its favorable Raven tracking conditions. Angles-only data 
obtained from another Raven telescope located at the AMOS Remote Maui 
Experiment (RME) facility was used to track the satellite of interest. Then the 
Analytical Graphics, Inc. Satellite Tool Kit Orbit Determination (STK/OD) program 
was used to compare the accuracy of the orbit prediction using ranging tracking 
data from Sirius and angles-only tracking data from Raven. This paper shows 
the improvement in orbit determination uncertainty obtained by adding Raven 
observations to the ranging data. The Raven angles data improved the orbit 
plane uncertainty and eccentricity estimate differences by over 80% when used 
with the range observations. 
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I. INTRODUCTION 


A. THE NEED FOR SPACE SURVEILLANCE 

The United States military’s Joint Publication 3-14 (Joint Doctrine for 
Space Operations) provides guidelines for planning and conducting joint space 
operations. The Joint Chiefs of Staff (JCS) recognize space as a significant force 
multiplier for the United States military due to the reliance on space systems to 
carry out operations. In addition, the civil sector’s reliance on space is a major 
factor in operations planning. Finally, it is recognized that current and future 
adversaries of the United States are dependent on their own space systems for 
intelligence collection against the United States. 

To that end, the Joint Publication 3-14 mandates that the United States 
establish space superiority: “The use of space control operations to support 
freedom of action in space will ensure the ability to provide space capabilities to 
the warfighter and deny the opposing force the same.” 1 In order to establish 
space superiority, four space missions must be accomplished: Space Control, 
Space Force Enhancement, Space Support, and Space Force Application. Since 
each of these mission areas is critical to achieving space superiority, they each 
have their own separate and distinct mission areas. Space Control operations 
“..provide freedom of action in space for friendly forces, while, when directed, 
denying it to an adversary, and include a broad aspect of protection of US and 
US allied space systems and degradation of adversary space systems.” 2 

Space Control encompasses the following mission areas: Protection, 
Prevention, Negation and Surveillance. The following (Fig. 1) from JP 3-14 
shows the relation of the four mission areas to space control 3 : 


1 Joint Pub 3-14, August 2002, p. IV-3. 

2 Ibid., p. IV-5. 

3 Ibid., p IV-7. 
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SPACE CONTROL MATRIX 


Gain or Maintain 
Space Control 



Battle Management Command, Control, 
Communications, Computers, and Intelligence 


Figure 1. JP 3-14 Space Control Matrix (From Ref. [1]) 


B. SPACE SURVEILLANCE 

While protection, prevention and negation are very important mission 
areas, this paper will focus on space surveillance. Space surveillance is deemed 
“fundamental to the ability to conduct the space control mission” and is defined 
as 


..requiring robust space surveillance for continual awareness of 
orbiting objects; real-time search and targeting-quality information; threat 

detection, identification, and location.conducted to detect, identify, 

assess, and track space objects and events to support space operations. 
Further, space situational awareness can be used to support terrestrially- 
based operations, such as reconnaissance avoidance and missile 
defense. 4 

It should be pointed out that space surveillance is critical to both providing 
freedom of action for friendly forces and denying freedom of action to enemy 
forces. It is the only mission area which affects both offensive and defensive 
sides of Space Control. By having a robust capability in space situational 
awareness, the United States will know when our forces are vulnerable to foreign 


4 Joint Pub 3-14, August 2002, p. IV-6. 
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intelligence-gathering space platforms in order to take timely and appropriate 
measures to defeat their attempts. 

C. SPACE SURVEILLANCE NETWORK 

As with any other surveillance technique or regime, the more accurate the 
data on the desired object, the better. In addition, minimizing the time and 
resources spent on obtaining highly accurate data is also a goal. This allows for 
a timely, clear picture of the battlespace, allowing the commanders more time to 
concentrate on avoiding or defeating the threat. 

Current space situational awareness systems (e.g., the Air Force Space 
Surveillance Network (SSN)) are meeting space situational awareness 
requirements, but are doing so at great cost and take thousands of personnel to 
operate and maintain the systems throughout the world. The Space Surveillance 
Network is comprised of over forty radar and optical sites throughout the world, 
with the majority located in the Northern Hemisphere. In 2001, the Air Force 
spent over $60 million to operate the Space Surveillance Network 5 . In today’s 
austere budget environment and with the age of the current SSN growing every 
year, it makes sense to investigate any and all ways to both increase accuracy 
and lower costs. 

D. HANDS 

One partial solution to this problem may be the High Accuracy Network 
Determination System (HANDS), a concept future network of optical telescopes 
that autonomously track both near-earth and deep space satellites and provide 
high accuracy orbit information. 6 Though HANDS is still in the development 
stage, the concept is to have thirty or more HANDS nodes spread throughout the 
world using automated telescopes to track earth orbiting objects of all regimes in 
conjunction with ranging data from selected SSN sites. 

5 Government Accounting Office Report GAO-02-403R, June 2002, p.2. 

6 Geosynchronous Orbit Determination Using HANDS, AAS 04-216, Sabol, Kelecy, Murai, 
February 2004, p. 1. 
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The angles data and ranging data are sent to the HANDS Operation 
Center by the Raven and AFSSN sites, where it is then fused and analyzed. The 
improved orbit estimates are then delivered to the customer. Fig. 2 below 
graphically illustrates the HANDS concept 7 : 



HANDS Ground 
Station (HGS) 

EAFB, RME 

•Satellite Imaging System 
•Site Command/Control 
•Weatherstation 


HANDS Operation 
Center (HOC) 

MHPCCSTACC 

•Data Fusion/Orbit Determination^ 
•Performance Metrics 
•RMA Metrics 
•Overall System Status 
•HANDS Command/Control 
•Space Situational Awareness 
•Space Object Identification 


Satellite Ranging 
Station (e.g. AFSCN) 

•Range Data Acquisition 


High Accuracy 
Orbit Determination 
SSA 
SOI 



Customer 


Figure 2. HANDS Architecture/CONOPS (From Ref. [3]) 


E. RAVEN TELESCOPE 

The RAVEN telescope is planned to be the backbone of the HANDS 
network. RAVEN is a class of small telescopes that combine inexpensive 
commercial hardware with state of the art astrometric image reduction 
techniques to produce high accuracy angular observations of satellites 8 . 

The RAVEN telescope system is comprised of five major components: the 
4 ft long telescope (0.37 meter mirror) and the dome which houses it, the 
telescope control computer, the Odin data processing workstation, the Global 
Positioning System (GPS) receiver and timing system, and a weather system 

7 Geosynchronous Orbit Determination using HANDS, AAS 04-216, Sabol, Kelecy, Murai, 
February 2004, p. 2. 

8 Ftecent Developments of the RAVEN Small Telescope Program, AAS 02-131, Sabol and 
others, January 2002, p.1. 
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which detects levels of wind, temperature, and humidity 9 . Raven telescope 
systems are designed to operate autonomously for weeks at a time without 
manual intervention 10 . 

This paper discusses the first trial of using the RAVEN telescope system 
for orbit determination accuracy against a HEO object. It will detail the object 
RAVEN tracked, the computer program used to determine Raven’s accuracy, 
and the results using methods to differentiate between a range only solution and 
a range plus angles solution. 


9 Ibid, p. 7. 

10 Recent Developments of the RAVEN Small Telescope Program, AAS 02-131, Sabol and 
others, January 2002, p.12. 
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II. BACKGROUND 


A. RAVEN TELESCOPE 

In the almost ten years since the Raven telescope prototype was built, 
Raven has become a very successful program. It has evolved from a good idea 
into a real-world Air Force Space Surveillance Network (SSN) sensor, capable of 
autonomously detecting, tracking and reporting geosynchronous objects with one 
arc-second accuracy. One Raven system is currently located at the summit of 
Mount Haleakala as a part of the Maui Space Surveillance Site (MSSS), 
providing space operators daily reports of deep space objects. A second system 
lies at the base of the mountain, in the Air Force Maui Optical and 
Supercomputing Site (AMOS) Remote Maui Experiment (RME) facility, hosting a 
large number of tracking experiments 11 . 

Due to Raven’s documented success in tracking geosynchronous (GEO) 
objects with sub arc-second accuracy 12 , it was decided to see how well it could 
track highly eccentric objects (HEO), which are typically more difficult to track 
and predict. FIEOs typically have an inclination of 63.4 degrees so that perigee 
does not drift, and usually perigee is in the southern hemisphere. Since the SSN 
radars are in the northern hemisphere the HEO satellites are usually beyond the 
detection range of the radars. Then their high latitude in the northern hemisphere 
makes it more difficult to obtain optical observations. Depending on the 
eccentricity the primary perturbations may be different at perigee and apogee. 
For those with a low perigee the high velocity at perigee means that atmospheric 
drag can have a significant effect as it passes through perigee. All of these 
factors combine to make their orbit determination and prediction more difficult. 
Generally speaking, if an object has an eccentricity greater than 0.1, it is 
considered highly eccentric. The Raven telescope’s accurate tracking of GEO 


^Recent Developments of the RAVEN Small Telescope Program, AAS 02-131, Sabol and 
others, January 2002, p.2. 

12 High Accuracy Orbit Analysis Test Results Using HANDS, Kelecy, Sabol, and Murai, Sept 
2003, p.9. 
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objects with one arc-second accuracy make it the perfect candidate to track the 
less predictable HEO objects. 


B. SIRIUS SATELLITE 

Though a Molniya-type orbit (e~0.7) is the most well-known of the highly 
eccentric orbits, access to current satellites flying in that orbit was not possible 
due to access to observations. To evaluate the effectiveness of the angle 
observations it is best to use a satellite with an accurate orbit. This allows one to 
compare the effectiveness of the orbit determination with angles only and with 
the addition of the angle observations to the primary observation set. The Sirius 
satellite radio constellation with an eccentricity of -0.27 satisfied this criterion. 
There are three satellites in the Sirius constellation. Sirius3, (Space Control 
Center number 26626) was chosen due to its long dwell time over the United 
States and favorable tracking conditions, that is, it stays illuminated by the sun 
while Raven is in umber. Sirius3 has a 24-hour period (16 hours in the northern 
hemisphere), and is commanded and tracked by a facility located in Quito, 
Ecuador (Lat 0.273 deg S, Long 281.5 deg E, altitude 2604 m). Below is a quick 
synopsis of Sirius3’s orbital elements and a ground trace (Fig. 3): 


Epoch: 09 Dec 2004 03:30:47.475 

Semi-major axis: 42165 km 
Eccentricity: .268 

True argument of latitude: 33.64 deg 
Inclination: 63.83 deg 

Right Ascension of the Ascending Node (RAAN): 29.82 deg 
Argument of Perigee: 269.76 deg 
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Educational Use Only 



Figure 3. Sirius3 Ground Trace (From STK 6.1) 


The 24-hour period and 63.4 degree inclination mean that the Sirius orbit 
is a double resonant orbit. 

Through Massachusetts Institute of Technology (MIT), nearly constant 
(over the time period vice every second of that time period) tracking and 
telemetry data for the Sirius3 satellite were obtained. Dr. Paul Cefola of MIT and 
Mr. Chris Croom from Sirius provided files with range and angle observations. 
These observations were then converted to B3 format, necessary for 
computation in the orbit determination software. There were approximately four 
months of Sirius observations to use. The observations included both range and 
angle (azimuth and elevation) portions and covered a time period from December 
4, 2004 through March 31, 2005. The Sirius angle observations are not very 
accurate, they help in the initial orbit determination, but their accuracy is not 
sufficient to improve in the maintenance of the orbit, i.e., the differential 
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correction. Consequently, they were only used in the initial orbit determination. 
Only ranging data were used in the analysis runs. 


C. STK/OD 

Analytical Graphic Incorporated’s (AGI) Satellite Tool Kit/Orbit 
Determination version 3.0 (STK/OD) was used for processing the observations 
and determining the orbit. The program was responsible for all of the “heavy 
lifting” of data processing. However, testing of the software was completed to 
ensure comparable results from separate approaches before serious attention 
could be used on the results. The STK/OD output was compared against the 
industry-standard Goddard Trajectory Determination System (GTDS) program’s 
output. Since GTDS is one of the recognized tools for orbit determination, 
settings in STK/OD were modified until output from STK/OD was similar to GTDS 
output. 


Many runs were made in trying to get the settings in STK/OD to provide a 
consistent result against GTDS. Orbital elements, covariance, residual plots, 
solar radiation pressure (SRP) plots (dynamic in STK/OD, fixed in GTDS), and 
position consistency plots were some of the products compared to ensure 
consistency between STK/OD and GTDS. Each run was built using the five main 
subsets of a STK/OD scenario: tracking facility, initial orbit determination, filter, 
smoother, and satellite. 

STK/OD requires that each scenario have the subsets listed above for 
every scenario run. The Initial Orbit Determination (IOD) is run first, using six 
angles (azimuth and elevation measurements at three different times) 
observations in order to compute a rough orbit. Once a rough orbit is generated, 
that orbit is transferred to the satellite. There is a least squares method that can 
be run to further refine the orbit with 10-20 of the initial observations, but it was 
decided that since the IOD was fairly close to the actual orbit, the least squares 
option did not need to be used. 


10 



1. Kalman Filter 

STK/OD’s filter is a forward-time recursive algorithm consisting of a 
repeating pattern of filter time updates of the state estimate, which propagates 
the state estimate forward, and filter measurement update of the state estimate 
which incorporates the next measurement. The filter uses the observations 
along with their location and a priori state estimate as the input, and provides 
optimal state estimates and realistic state error covariance matrices as the 
output, updated after every observation and at 1-second intervals 13 . The initial 
covariance is input as “orbit uncertainty”, listed in the satellites settings. It is 
input in the radial, in-track, cross tack (RIC) reference frame. The only 
requirement for the initial covariance is that it not be too small, with many 
observations the final covariance is essentially independent of the initial 
covariance as long as it is not too small. Typically, it should be at least an order 
of magnitude larger than the expected final covariance. The only effect of a larger 
initial covariance is that the time to converge to a “steady state” covariance 
increases. For the runs, a diagonal covariance was used with 100,000 m for the 
RIC standard deviations and 100 m/sec for the RIC rate standard deviations. 

2. Process Noise 

In a Kalman filter implementation, process noise is used to represent the 
unmodeled accelerations and to prevent the covariances from getting too small. 
STK/OD has three types of process noise. To capture the gravitational force 
uncertainties and other unknown forces such as outgassing there is process 
noise in the radial, in-track and cross track directions (RIC). The magnitude in 
each of these three directions is an input quantity. In addition, there is process 
noise associated with both the solar radiation and atmospheric drag forces. 
Since the orbit’s perigee was well above the atmospheric drag region (-29,500 
km), atmospheric drag was not considered. The key question is what should the 
standard deviation of the process noise be, particularly with the different types of 
process noise. Due to the fact that a high order gravity model was used and lunar 
13 STK/OD manual. 
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and solar perturbations were included it was expected that the solar radiation 
pressure process noise would capture all the uncertainty. For the solar radiation 
the satellite was modeled as a sphere. To allow for the fact that the satellite is 
not a sphere an uncertainty in the two directions perpendicular to the satellite-sun 
line is allowed. This uncertainty is modeled as process noise with a magnitude of 
0.3 times the solar radiation acceleration along the satellite-sun line with a half- 
life of 300 minutes. Including only solar radiation uncertainty was not sufficient, 
the residuals were too large. This was not due to station keeping maneuvers as 
a maneuver schedule was provided with the data and no maneuver was 
performed during the analysis time. Possibly some outgassing was occurring or 
there were momentum dumps that resulted in small linear accelerations that are 
caused by thruster mismatch. Therefore, it was necessary to include process 
noise in the RIC directions. A balance needed to be found between too much 
and too little unmodeled process noise in the satellite settings, and through 
numerous trials, it was found that a process noise of 0.01 cm/s for each of the 
three RIC directions resulted in the necessary consistency in the residuals. 
Consequently, a process noise value of 0.01 cm/s was used for each axis. 
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III. ANALYSIS 


A. OBSERVATIONS 

A discussion of the observations is necessary before the analysis is 
presented. As discussed previously, observations were provided to AMOS from 
Sirius satellite radio via MIT. These observations were converted to the B3 
format required for STK/OD. The Sirius observations in the month of December 
were fairly consistent, with an average of six observations per hour, spaced 
within a two minute timeframe, usually between :30 and :32 minutes past the 
hour, Greenwich Mean Time (GMT). An example of the raw Sirius observation 
data is below: 


Source Ant Type Status Data 

2004/12/09 15:30:22.264 c3_120920041530 QTB range accept 39325.0201 km 

2004/12/09 15:30:24.300 c3_120920041530 QTB azimuth accept 305.698000 degs 

2004/12/09 15:30:24.300 c3_120920041530 QTB elevation accept 43.8500000 degs 


Estimate 
39367.1028 km 
305.059170 degs 
43.5886859 degs 


Noise 

10.0000000 meters 
0.0200000000 degs 
0.0200000000 degs 


Residual 
-42082.6177 m 
0.638830439 degs 
0.261314137 degs 


The December Sirius observations were fairly consistent, except for some 
timeframes which were missing observations. Unfortunately, some of these 
timeframes coincided with the timing of the Raven observations. Fig. 4 below 
illustrates a timeframe in which there were gaps in Sirius observations (indicated 
by the TrackerlD 999.00), coincident with Raven observations (shown as 998.1 
and 998.2): 
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Measurement Times 


6750 
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8750 


9750 
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17750 18750 19750 20750 21750 


Range 


Figure 4. Sirius and RME Raven Observation Times 

Though there were many observations collected by the AMOS Raven, 
located at the summit of Mount Haleakala, from December into January, these 
observations were not useable, due to shutter control issues. That left 
observations from the RME Raven, which was able to get approximately 200 
observations of azimuth and elevation during the 15-17 December timeframe. 
These observations are the basis of the analysis 14 , and there were enough 
observations to give a fairly clear picture of the results. 


The majority of RME Raven observations were collected on 15-16 
December, with ten observations of azimuth and elevation collected on 17 
December. However, due to the exceptional one arc-second accuracy of the 
Raven telescope, the lack of observations is compensated by their accuracy. 


14 Observations made at the Maui Space Surveillance System (MSSS), Maui, Hawaii, USA 
are the result of collaboration between Dr. Chris Sabol and Detachment 15 of the US Air Force 
Research Laboratory, which owns and operates the MSSS. 
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B. DATA QUALITY TESTS 

Once the observations were collected, an orbit had to be established in 
order to have a way of measuring improvement while determining the correct 
settings in STK/OD. The baseline orbit was chosen to be the Sirius orbit from 9- 
17 December, because of the number and quality of the range observations from 
the satellite. A number of tests were performed on the data to ensure the 
settings in STK/OD and the quality of observations were sufficient to establish 
the baseline orbit. One of these tests was the residuals check, shown in the 
graph in Fig. 5. A sigma of 10 meters was used for the Quito tracker. The 
residuals were consistent and reflected the expected result: 
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Figure 5. Sirius Range Only Residuals 


Another test that was completed was to examine the Solar Radiation 
Pressure (SRP) dynamics plot. In STK/OD, the SRP estimate is shown as a 
correction to the estimated (input) value and the filter estimate. Since the Sirius 
orbit has an apogee of -54,800 km, solar pressure on the satellite is a factor in 
orbit determination. To ensure an accurate solution, SRP was examined for any 
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perturbations which may skew the results. The figure below shows that SRP was 
consistent enough as to not affect the results. Both the SRP estimate and the +/- 
2 sigma values are shown in Fig. 6. 



Figure 6. Solar Radiation Pressure Estimate 


A final test for goodness was the comparison of the filtered run versus the 
smoothed run. The smoothed run is obtained by using the state and covariance 
at the final observation and performing a “backwards” filter with all the 
observations to obtain the optimal estimate during the span of data. This position 
consistency test ensures that the filter behaves as expected and there was not a 
large difference in the filtered results and the smoothed results. Fig. 7 shows the 
normalized differences of range, in-track, and cross-track results. The RIC 
values depicted in Fig. 7 show the difference in the filter run and the smoothed 
run using the McReynolds consistency test 15 found in STK/OD. The mean is 
zero, while the upper and lower bounds are +/- 3 sigma (dimensionless), visible 
at the upper and lower edges of the graph. Values within the +/- 3 sigma are 
considered reasonable. 

15 Optimal Orbit Determination, STK/OD White Paper, Wright, J., p.5. 
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Figure 7. Position Consistency Graph 


After establishing this orbit as the baseline orbit, the addition of the Raven 
observations and comparing the outputs of each was the next step. All settings 
in STK/OD were kept the same as they were for the Sirius observations only 
runs, with the exception of settings which added the Raven azimuth and 
elevation observations to the scenario. These settings can be seen in Appendix 
A. 


1. Orbit Determination Methods 

A measure of the quality of an orbit determination is the accuracy of the 
orbit prediction. The two primary approaches for determining this quality for 
different scenarios and sets of measurements are: (1) comparison of the actual 
orbit determination results, and (2) comparison of the covariances. 

Comparing the accuracy of actual orbit determination results for different 
scenarios and sets of measurements of orbit determination is difficult when there 
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is no concrete truth orbit. One approach often used is to perform an “abutment 
check”, perform an orbit determination for two separate fit spans, then propagate 
the state from one epoch of one fit span across the other fit span, or if the fit 
spans are not contiguous propagate each orbit to a common time point between 
the two fit spans, and difference the results. This is shown graphically in Fig. 8 
below: 


Fit Span 1 

Fit Span 2 

III! iftibiili 

ig iiiillllllliillllliillll 

Orbit 1 

: 

■ 

Orbit 1 

Orbi 

. . . . 

88888 

t Compar 

ison 


Figure 8. Abutment Check 


If the orbit for the second fit span is truth, the difference is the error. If it is 
not truth, the problem becomes determining how to interpret the results. In the 
in-track direction, there is secular growth so if you predict long enough, this error 
will dominate the other errors during the fit span and one can assume the orbit 
during the fit span is truth. This is not the case for radial and cross-track 
comparisons because these errors are periodic and generally do not grow. 
Consequently, the comparison of real world results usually only works in the in¬ 
track direction. Since the “abutment check” is a comparison of actual results each 
case represents only one sample. To make any real conclusion from the 
“abutment check” there needs to be enough cases for a statistically reliable 
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sample. Unfortunately, the lack of angle observations over an extended period of 
time resulted in only one case, which prevented any meaningful comparison of 
the actual results in the in-track direction. Thus, the approach focused on 
comparing elements of the covariance. 

2. Covariance 

The covariance provides information on the accuracy of orbit 
determination. It can be used for comparison or just for a single fit. Of course, 
for the covariance to provide a valid assessment of the orbit determination 
accuracy, the modeling and sensors errors have to be accurately modeled. If 
one is comparing covariances, one can reasonably argue in some cases if both 
are in error in the same manner, the comparison is valid. For example, if the 
actual measurement errors for both fit spans are 5 arc-seconds, but are modeled 
as 10 arc-seconds, then the comparison should be valid. 

After the last time step in the filter run, STK/OD outputs a covariance 
matrix in the radial, in-track, and cross-track (RIC) reference frame. The 
following are the filtered covariance and orbital elements outputs from the Sirius 
range only, and from the Sirius range plus Raven angles observations from the 
December 9-17 fit span: 


Sirius Ranqe-Onlv: 


Final Value 


Semimajor axis (km) 


42164.509729 

Eccentricity 


0.268586659 


True Arg of latitude (deg) 


350.4501624 


Inclination (deg) 


63.8292837 


Right Ascension of Ascending Node (RAAN) (deg) 29.6918866 


Arg of Perigee (deg) 


269.8341598 


RIC Sigma Correlation Matrix (m & cm/s) 



83.94 m -0.86 0.92 

0.89 

-0.91 

0.88 

660.34 m -0.95 

-1.00 

0.98 

-0.92 

546.94 m 

0.95 

-0.98 

0.96 


5.90 cm/s 

-0.98 

0.92 



2.15 cm/s 

-0.95 

2.22 cm/s 
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Here the following matrix is represented: 


PxY 

P X z 

P XX 

P AT 

P XZ 

(J Y 

P Y Z 

Pyz 

P YY 

Pyz 



P zx 

P ZY 

P zz 



°x 

P AT 

P XZ 





Pyz 





°z 


Here, a x , o Y , and a z denote the standard deviations of the radial error 
(X), in-track error (Y), and cross-track error (Z) and their rates of change 
[a t ,(7 Y ,(7 z ). p values are the correlation coefficients of the two values listed in 
the subscript of each p. 

Making direct comparisons of these quantities provides no useful 
information because they are periodic and possibly secular. The in-track error 
growth rate (drift rate) is caused by a semi-major axis error. The radial error is 
caused primarily by an eccentricity error and the cross-track error is caused by 
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an error in the orbit plane estimate. The following equations 16 are taken from 
reference (7). The drift rate, d, is proportional to the semi-major axis error: 


( 1 ) 


where da is the semi-major axis error. The orbit plane error /is given by 


y 2 = Si 2 + (3Q.sm i ) (2) 

where (i,Si) are the inclination and inclination error, and, is the right 

ascension error. From ref [7], these quantities as a function of the RIC errors and 
error rates are 


where 


da = — f (esinf)X + (l + ecos f)Y ~] + —- ——^—X 

nr/ L J rf 


(2tJ 4 ] 

f fa > 

--1 
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(y, 2 +yj) 
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2V„ 
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( 3 ) 


p = a( 1 - e 2 ) 

R — p/( 1 + ecos/) 

V = (l-e 2 ) 1/2 

A = Jup (4) 

^ =(l + ecos f)(J(iuTp)) 

V r = (esin f)(^J(p / p)) 

V c =JjuJa 


16 The State Transition Matrix of Relative Motion for the Perturbed Non-Circular 
Reference Orbit, by Gim, D. and Alfriend, K., AIAA J. of Guidance, Control, and Dynamics, Vol. 
26, No. 6 November-December 2003, pp 956-971. 
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In equations (3) and (4), {X,X) = radial, radial rate, ( Y,Y) = in-track, in¬ 
track rate, (Z,Z) = cross-track, cross-track rate. V t ,V r ,V c are the tangential, radial 
and circular velocities, respectively, p is the semi-latus rectum, R is the radius, 
a is the semi-major axis, h is the angular momentum, n is the gravitational 
parameter, and / is the true anomaly. 

The standard deviations of the quantities in equation (3) are 


Semi-major Axis Error: 



v) 2 M ! 


a a 


— + — — + - — + 

K [V [V (r ( R 


K Y <r v Y a 


K K V 


a I G x I Vy) o 
- — — Ay+2 

RAvIrP* 


K Y a] a 


KIR V 



( 5 ) 
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Orbit Plane Error: 


°r = 
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C. RESULTS 

Using the STK/OD provided classical orbital elements and covariance 
matrices, the numerical values for semi-major axis error, orbit plane error, and 
eccentricity error for the range only and the range and angles cases are 
summarized in Table (1) below. 


Metric 

Semi-maior-Axis Error (m) 
Orbit Plane Error (deg) 
Eccentricity Error 


Range Only 
109.226 
0.000857 
0.0000283 


Ranqe+Anqles 

108.754 

0.000127 

0.00000539 


Difference Imp % 

0.472 0 . 432 % 

0.000730 85 . 221 % 

0.0000229 81 . 028 % 


Table 1. Metric Results Comparison 


What can be seen here is that there is a significant improvement of over 
80% in both the orbit plane uncertainty and eccentricity estimate differences 
when Raven angles observations are used in conjunction with the range 
observations, even though there is only a very slight improvement in the semi¬ 
major axis error. This demonstrates that Raven angle observations improved the 
orbit determination parameters for this satellite. 
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IV. CONCLUSIONS 


The Raven telescope is a valuable asset to the space surveillance 
network, as it has already proved itself in tracking Geosynchronous satellites. 
Due to its accuracy, it was a natural fit to attempt to track (with the same orbit 
determination accuracy) Highly Eccentric Orbiting satellites. This study, though 
limited, has shown that Raven is ready to take the next step in completing the 
goal to track any deep space man-made object orbiting the earth. By 
strengthening the quality of the HEO object’s orbit parameters, follow-on tracking 
is improved with a smaller search area for the object. This leads to better space 
situational awareness. 

Because there was only three days of angles data, it was not possible to 
completely verify Raven could track HEO objects with the same accuracy. For 
this thesis, it would have been better to have a full week or more of Raven angle 
observations in which to build an orbit from, then to compare with a range only 
orbit. Abutment checks could have been completed and used for another proof 
of the increased orbit determination accuracy of the Raven telescope. However, 
based on the improvement in critical orbit metrics, this analysis makes a strong 
point for follow on testing against more objects in highly eccentric orbits. 

For a complete evaluation of the angle observation contribution, a solid 
30-day period of Raven angles observations would be enough to not only add to 
range only observations as done in this thesis, but to build an entirely new orbit. 
This new angle-only orbit could then be compared to a range-only orbit. If the 
two orbits compare favorably, it would indicate that the Raven could be a reliable 
sensor to track Highly Eccentric Orbits. There is more work to be done and all 
indications are that with enough data and sufficient documentation, HEO objects 
could then be added to Raven’s object tracking list. 
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APPENDIX A. STK/OD SCENARIO SETTINGS 


1. FACILITY (Quito) 

• Position Geodetic 

o Lat -0.273042 deg 

o Lon 281.524 deg 

o Alt 2604 m 

• Tracking ID 999 

• Estimate Nothing 

• MinElevation 5 deg 

• MaxElevation 90 deg 

• RangingMethod Transponder 

• AntennaType Mechanical 

• Optical Properties 

o PolarExclusion 1 deg 

o ReferenceFrame MEME of Date 

o AberrationCorrections None 

• TroposphereModel 

o Enabled No 

o Model SCF 

• TroposphereData 

o SurfaceRefractivity Constant 

o Value 340 

• lonosphereModel 

o Enabled No 

o Model IRI2001 

o TransmitFreq 2267.5MHz 

o ReceiveFreq 1815.77MHz 

2. FACILITY (RMERaven) 

• Position Geodetic 

o Lat 20.7462 deg 

o Lon 203.568 deg 

o Alt 105.38 m 

• Tracking ID 998 

• Estimate Nothing 

• MinElevation 5 deg 

• MaxElevation 90 deg 

• RangingMethod SkinTrack 

• AntennaType Optical 

• Optical Properties 

o PolarExclusion 1 deg 

o ReferenceFrame MEME J2000 

o AberrationCorrections None 

• TroposphereModel 

o Enabled No 

o Model SCF 

• TroposphereData 

o SurfaceRefractivity Constant 

o Value 340 

• lonosphereModel 

o Enabled No 

o Model IRI2001 

o TransmitFreq 2267.5MHz 

o ReceiveFreq 1815.77MHz 
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3. INITIAL ORBIT DETERMINATION (IOD) 


• Method GoodingAnglesOnly 

o TrackerList Quito 

o StartTime 09 Dec 2004 00:00:00.000 UTCG 

o StopMode LastMeasurement 

o MeasurementSampleSize 300 

o MinimumElevation 5 deg 

o SelectedMeasurments Double click to edit 
o HalfRevEstimate 0 

o Lambertlndicator 0 

o Rangel Estimate 5 Re 

o Range3Estimate 5 Re 

o Maxlterations 25 

o ConvergenceValue 1e-012 

o HalleyNewtonLimit0.5 

o NumericPartialEpsilon 1e-005 

o T12 7176.33 sec 

o T13 10742.2 sec 

• Solutions 

o NumberOfSolutions 2 

o UseSolution 1 

• Output 

o OrbitState Keplerian 

o CoordinateFrame J2000 


4. FILTER 

• ProcessControl 

o StartMode Initial 

o StartTime 9 Dec 2004 03:30:47.475 UTCG 

o StopMode StopTime 

o StopTime 31 Dec 2004 23:59:59.990 UTCG 

o ProcessNoiseUpdatelnterval 1 min 

• Restart 

o SaveRecordstoFile false 

o MaxRecordsinFile 100 

o SaveFrequency 60 min 

• OptionalSolveForParms 

o MeasBiases true 

• Output 

o DataArchive 

■ OutputStateHistory AIITimes 

■ EveryNSteps 1 

■ SaveOnlyLastMeasPerStep false 

■ OutputMeasHistory true 

■ OutputManeuvers false 

■ OutputSummary true 

■ OutputHistograms true 

■ HistogramSize 3 

• NumberHistorgramBins 22 

o Display 

• EveryNMeasUpdates 1 

■ EveryNTimeUpdates 1 

■ ShowPassTimes true 

o SmootherData 

■ Generate true 

■ TimeMode FilterSpan 

o STKEphemeris 

■ DuringProcess 

• Generate false 

■ Predict 
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Generate 


false 


5. SMOOTHER 

• Input 

o Files double click to edit 

o Remove false 

• ProcessControl 

o StartMode LatestFilterTime 

o StartTime 31 Dec 2004 23:59:59.990 UTCG 

o StopMode EarliestFilterTime 

o StopTime 9 Dec 2004 03:30:47.475 UTCG 

o OutputLag 0 min 

o interval Length 1440 min 

o IntervalOverlap 720 min 

• Output 

o DataArchive 

■ OutputStateHistory AIITimes 

■ EveryNSteps 1 

■ OutputManeuvers true 

o STKEphemeris 

■ DuringProcess 

• Generate true 

• TimeGrid Filter 

■ Predict 

• Generate true 

• TimeStep 1 min 

• StopMode TimeSpan 

• TimeSpan 720 min 

■ Covariance true 

■ CovarianceType position 3x3 Covariance 


6. SATELLITE 

• Description 

• OrbitState Keplerian 

• EstimateOrbit true 

• OrbitClass LOeHEO 

• PhysicalProperties 

o Mass 3000 kg 

• Measurementprocessing 

o TrackingID 26626 

o MeasurementTypes Range Azimuth Elevation 

o ResidualEditing 

■ NominalSigma 3 

■ Dynamic 

• Enabled true 

• HighSigma 10 

• NumRejectToStart 2 

• NumAcceptToStop 10 

• InitialHighSigmaDuration 120 min 

o ThinningTime 0 sec 

o MinPassDelta 20 min 

• MeasurementStatistics None 

• MinGrazingAlt 100000 m 

• OpticalProperties 

o PolarExclusion 1 deg 

o ReferenceFrame MEME of Date 

o AberrationCorrections None 

• RangingMethod Transponder 

• lonosphereModel 

o Enabled false 

• ForceModel 

o Gravity 

■ DegreeandOrder 12 
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Tides 


o 


o 


• SolidTides false 

• OceanTides false 

■ GeneralRelativityCorrection false 

■ VariationalEquations 

• Degree 2 

■ ProcessNoise 

• Use BasedOnOrbitClass 

• WillUseProcessNoise true 

• OmissionErrorModeling 

o Enabled false 

o Scale 1 

• CommissionErrorModeling 

o Enabled false 

o Scale 1 

• ThirdBodies 

o Sun true 

o Moon true 

o Planets flase 

o UseinVariationaEquations false 

Drag 


■ Use BasedOnOrbit 

■ WillUseAirDrag false 
SolarPressure 

■ Use BasedOnOrbit 

■ WillUseSolarPressure 

■ EstimateSRP 

■ CPNominal 

■ Area 

• CPInitialEstimate 

■ CPHalfLife 

■ ReflectionModel 

■ SunPosMethod 

■ UselnVariationalEquations 

■ AddProcessNoise 

■ EclipticNorthFraction 

■ EclipticPIaneFraction 


true 

true 

3 

38 m A 2 
0 

300 min 

Sphere with diffuse reflection 

ApparentT oT rueCB 

true 

true 

0.3 

0.3 


o Plugin 

■ Use false 
o UnmodeledAccelerations 

■ ProcessNoise 

• RadialVelocitySigma 

• IntrackVelocitySigma 

• CrosstrackVelocitySigma 

• Timelnterval 

■ InstantManeuvers 

■ FiniteManeuvers 

■ OrbitErrorTransitionMethod 
PropagatorControls 

o IntegrationMethod RKF 7(8) 

o StepSize 


0.01 cm*sec A -1 
0.01 cm*sec A -1 
0.01 cm*sec A -1 
2 min 

InstantManeuvers 

FiniteManeuvers 

VariationalEquations 


■ Time .5 min 

■ TrueAnomaly 2 deg 

■ EccentricityThreshold 
EphemerisGeneration 

o CreateSTKFile false 

OrbitUncertainty 
o RSigma 


0.04 


o ISigma 
o C_Sigma 
o Rdot_Sigma 
o ldot_Sigma 
o Cdot_Sigma 
o AllCorrelations 
FilterEvents 

o MeasurementRejectThreshold 

■ NumForWarning 0 

■ NumForAlert 0 

o MeasurementAcceptTimer 


100000 m 
100000 m 
100000 m 
100 m*sec A -1 
100 m*sec A -1 
100 m*sec A -1 
0 
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TimeGapForAlert 


TimeGapForWarning 

0 


0 min 
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